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Abstract: We study, globaly in time, the velocity distribution f(v,t) of a spa- 
tially homogeneous system that models a system of electrons in a weakly ionized 
plasma, subjected to a constant external electric field E. The density / satisfies 
a Boltzmann type kinetic equation containing a full nonlinear electron-electron 
collision term as well as linear terms representing collisions with reservoir parti- 
cles having a specified Maxwellian distribution. We show that when the constant 
in front of the nonlinear collision kernel, thought of as a scaling parameter, is 
sufficiently strong, then the L 1 distance between / and a certain time dependent 
Maxwellian stays small uniformly in t. Moreover, the mean and variance of this 
time dependent Maxwellian satisfy a coupled set of nonlinear ODE's that consti- 
tute the "hydro dynamical" equations for this kinetic system. This remain true 
even when these ODE's have non-unique equilibria, thus proving the existence of 
multiple stabe stationary solutions for the full kinetic model. Our approach relies 
on scale independent estimates for the kinetic equation, and entropy production 
estimates. The novel aspects of this approach may be useful in other problems 
concerning the relation between the kinetic and hydrodynamic scales globably in 
time. 



1. Introduction 

The mathematical understanding of equilibrium phenomena has greatly advanced in 
the past few decades. The elegant and precise theory of Gibbs measures provides a di- 
rect bridge between the microscopic and macroscopic descriptions of such systems. This 
includes a general conceptual framework as well as nontrivial explicit examples of the 
coexistence of multiple equilibrium phases for certain values of the macroscopic control 
parameters. 

There is no comparable general theory for nonequilibrium systems and the micro- 
scopic study of phase transition phenomena in such situations appears to be far beyond 
our mathematical grasp at the present time. Our mathematical understanding of the great 
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variety of nonequilibrium phase transitions observed in fluids, plasmas, lasers, etc., there- 
fore depends entirely on the study of bifurcations and other singular phenomena occurring 
in the nonlinear equations describing the macroscopic time evolution of such systems. 

While there has been much progress recently in deriving such equations from simple 
microscopic and even realistic mesoscopic model evolutions the passage to the macroscopic 
scale is well understood only over time intervals in which the solutions of the macroscopic 
equations stay smooth. This is true for example in the passage from kinetic theory, where 
the evolution is described by the Boltzmann equation, to hydrodynamics, where it is 
described by either the compressible Euler or Navier-Stokes equations, depending on how 
we choose our macroscopic time scale [1, 2, 3, 4, 5, 6]. These derivations, which are 
based on Chapman-Enskog type expansions, require for their validity the uniqueness and 
smoothness of the solutions of the hydrodynamic equations. The reason for this is that 
control of the error terms in the estimates depends on a-priori smoothness estimates for 
solutions of the macroscopic equations. Thus, they shed no light on the actual behavior 
of the mesoscopic description when the solution of the hydrodynamical equations develop 
singular behavior. 

To overcome this problem it is clearly desirable to develop methods in which one does 
not use any a-priori smoothness estimates for solutions of the macroscopic equations, but 
instead uses scale independent estimates on the mesoscopic equation. This is what we do 
here for a simple model inspired by plasma physics [7, 8, 9]. 

Our starting point is a kinetic theory description of the system. Grave difficulties are 
posed by the fact that as of yet, not very much is known in the way of a-priori regularity 
estimates for solutions of the spatially inhomogeneous Boltzmann equation. This is quite 
different however, from the lack of estimates for the macroscopic equations - there it is 
clear that in the interesting cases the desired estimates just don't exist. Shock waves do 
form. In the Boltzmann case however, it is likely that a-priori regularity estimates in the 
velocity variables, say, invariant under the Euler scaling, are there, but simply have not 
yet been discovered. Still, the lack of such estimates is a grave difficulty in the way of 
rigorous investigation of the problem at hand. 

We sidestep this difficulty by considering a spatially homogeneous system, but one that 
is driven by an electric field, and coupled to heat reservoirs. In this case the usual hydro- 
dynamic moments are not conserved and the system will have non-equilibrium stationary 
states. We prove then in a certain simplified, but still recognizable physical situations, that 
the kinetic description closely tracks the macroscopic description even when the driving 
is sufficiently strong for the latter to undergo phase transitions. More precisely, we show 
that the velocity distribution function is close to a Maxwellian parametrized by a tem- 
perature T and mean velocity u which satisfy certain non-linear equations, which are the 
macroscopic equations for this system. Moreover, it does so globally in time, even when 
the stationary solutions of these macroscopic equations are nonunique. 

We are in fact particularly concerned with the stability of of these stationary solutions 
- the existence of multiple stationary states being analogous to the coexistence of phases 
in equilibrium systems. For such questions we need results that guarantee that a solution 
of the kinetic equations will stay near a solution of the macroscopic equations globaly in 
time. This seems to be difficult to accomplish by standard expansion methods, at least 
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in the range of driving field strengths where the macroscopic equations have the most 
interesting behavior. Instead of expansion methods, we use entropy production [10, 11] 
to show that the solution of the kinetic equations must stay close to some Maxwellian, 
globally in time. Then, we show that the moments of this Maxwellian must nearly satisfy 
the macroscopic equations. In this way we get our results. The next section specifies the 
model more closely, and states our main results. A preliminary account of this work in 
which the Boltzmann collisions were modeled by a BGK collision kernel was presented in 
[12]. 

2. The model and the results 

Our formal set up is as follows: We consider a weakly ionized gas in IR 3 in the 
presence of an externally imposed constant electric field E. The density of the gas, the 
degree of ionization and the strength of the field are assumed to be such that: (i) the 
interactions between the electrons can be described by some nonlinear, Boltzmann type 
collision operator, and (ii) collisions between the electrons and the heavy components of 
the plasma, ions and neutrals, are adequately described by assuming the latter ones to have 
a spatially homogeneous time independent Maxwellian distribution with an a proiri given 
temperature [7]. Under these conditions the time evolution of the spatially homogeneous 
velocity distribution function f(v,t) will satisfy a Boltzmann type equation 

= —E ■ V/ + Lf + e" 1 Q(/), (2.1) 

where V is the gradient with respect to v in IR 3 , E is a constant force field and Q is a 
nonlinear collision term which will take either the form of the Boltzmann collision kernel 
for Maxwellian molecules, or the one corresponding to the BGK model. We treat both 
cases here because it is possible to provide a little more detail concerning the nature of 
the equilibria in the BGK case. The parameter e > is thought of as a scaling parameter 
that goes to zero in the hydrodynamical limit. The linear operator L represents the effect 
of collisions with reservoir particles. It will be assumed to have the form: 

Lf{v)=L 1 f{v)+L 2 f{v), (2.2) 

with 

Lif(v) = V • (D(v)M(v)v(j^y), (2.3) 

a Fokker-Planck operator, representing energy exchanges with the reservoir assumed to be 
at temperature T = 1, so that 

M(v) = (27r)- 3/2 exp(-H 2 /2) 

and 

D(y) = aexp(-%| 2 /2) + c (2.4) 
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for some strictly positive constants a, b and c: the symbols a, b, c shall henceforth always 
refer to these parameters wherever they appear. The specific form (2.4) of the velocity 
space diffusion coefficient is not important. We specify it for sake of concreteness. The 
properties we really need for D(v) will be clear from the proofs. The operator L2 represents 
momentum exchanges with the heavy reservoir particles and is given by 

L 2 f(v) = v[f{v) - /(«)], (2.5) 

with v a positive constant and f(v) the sphericalized average of f(v). 

For any probability density /, we shall let Mf denote the Maxwellian density with 
the same first and second moments as /. Explicitly, 



M / = (2 7 rT)- 3 / 2 exp[-^^], 



with 



(2.6) 



u:= / vf(v)d 6 v, (2.7) 

'IR 3 



U v2 k 

1 Jtr 3 



(v)d 6 v, (2.8) 

z Jtr 3 

and T = |(e — \u 2 ). In the BGK model the collision kernel is 

QBGK(f)=M f -f . 

The Boltzmann collision kernel is given by 

Qb(/)(v)= f dvJ duB{\v - v*\,u) [iV)/«) - f(v)f(v*)] . (2.9) 
Jtr 3 Js+ 

Here S£ = {u G R 3 : u 2 = 1, u ■ (v - u*) > 0}, and 



v =v — u) • (v — v*)u, 



(2.10) 



are the outgoing velocities in a collision with incoming velocities v and v* and impact 
parameter uj, B(\v — v#\,u) is the collision cross section, depending on the intermolecular 
interactions. For Maxwellian molecules, with a Grad angular cut-off [Gr], 

B(\v-v*\,u) = h(6), (2.11) 

with 9 the azimuthal angle of the spherical coordinates in £2 with polar axis along v — i>* 
and h(9) a smooth non negative bounded function. Thus, for any normalized / we can 
write 

Q B (f)=£(fof-f) (2.12) 
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with 



f°f( v )=hf dv * f d9h(9)\ S m9\f(v')f(v'^ 
liL Je 3 Jo Jo 



£=tt [ deh(6)\sm6\ > 0, 
Jo 



(2-13) 



With either Qbgk or Qb for the collision kernel in (2.1), this term tends to keep / close 
to Mf, and one could certainly expect this effect to dominate for small values of e. Thus, 
formally in the limit as e vanishes, / will actually equal Mf for all time t, and to keep 
track of its evolution, we need only keep track of u(i) and e(t). 

Using the prescription / = Mf, in the right side of (2.1), one easily evaluates the time 
derivative of the first two moments of the so modified (2.1), to obtain formally 



d(u(t)\_ 
dt I e(t) 



F(u(t),e(t)) 
G(u(t),e(t)) 



(2.14) 



The functions F and G are given explicitly by 

F(u, e) = E — u[u + c + a exp(— w) 
aexp(—w 



1 + b 



(1 + 6//?) 5/2 



(2.15) 



G(u, e)=Eu-c 2e(l - 0) + 0u 



(1 + b/ [If* 



1-P, 



2e ( 1 -(3) +u ^(3- b —J-) , (2.16) 



where (3 = T _1 and w = bu 2 /2(l + bT). The tildes in (2.14) are to remind us that the 
equation is valid only when / = Mf. 

Equations (2.14) represent the hydrodynamical description of the gas. Our primary 
goal here is to show that such a description actually does hold for small, but positive, 
values of e, i.e. that the interaction between the hydrodynamic and the non-hydrodynamic 
modes does not destroy the picture involving only the hydrodynamic modes. The following 
theorem enables us to do this. 



Theorem 2.1. Let f be a solution of (2.1) with 



/(.,0) = M /( ., ) • 

Then, for any fixed integer k > 0, there is an €q > and functions 5i(e), 62(e), going to 

zero as e — > 0, depending only on a, b, c, \E\ and e(0), such that for e < eo the solution of 
(2.1) satisfies 

sup \\f(-,t) - M /( ., t) \\ LHM s } < 8 1 (e) (2.17) 

and 

sup || V fc (/(-,t)-M /( ., t) ) 11^3) <5 2 (e) (2.18) 

for k < ho- 
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Remark: The assumption on /(•, 0) is not essential: our methods allow us to easily modify 
the result to take into account an initial layer. 



We shall use Theorem 2.1, in the proof of Theorem 2.2 below, that if we compute 
u(t) and e(t) for a solution / of (2.1) satisfying the conditions of Theorem 2.1, then the 
moments of / will satisfy the equation, 

d ( U (t)\ _ (F(u(t),e(t))\ 1/2 / 7 (f)\ 

dt {e(t)J ~ {G(u(t),e(t))) +d{€ > { V (t)J {2 - W > 

where F and G are the non-linear functions of u and e given by (2.15) and (2.16), that arise 
in the e = limit, and 7 (t) and r](t) are bounded uniformly in t with a bound independent 
of e. Of course they depend on the full solution / of (2.1), but such estimates and a simple 
comparison argument will then lead to the following theorem, which says that the system 
(2.14) obtained in the e — > limit does give an accurate picture of the small e regime. 

Theorem 2.2 Let (-u*,e*) be a stable fixed point of the system (2.14) and let M( u * e *) be 
the corresponding Maxwellian density, with moments u* and e* . Then given any 5 > 0, 
there is an e greater than zero such that i// e (-,t) solves (2.1) with this value of e and 

||/ e (.,0)-M (u . ie . ) (.)|| L i( R 3) <6 (2.20) 

then 

||/ e M)-M (u . ie .)(.)||Li(]R3) <<5 (2-21) 

for all t > 0. 

//, however, (w*,e*) is not stable, then there exist a 5 > so that for every e > 0, 
there is a solution / e (-,t) of (2.1) with Maxwellian initial data satisfying 

||/ e (.,0)-M (li . )e . ) (.)|| L i(R3) <e, (2.22) 

but such that for some finite t > 

||/ e (-,*)-M (tt . ie , ) (.)||Li(R3 ) >S . (2.23) 



The proof of these theorems, which are fairly complicated even for the simple BGK 
model will be given in the next sections. 

The above theorems allow us to rigorously prove that our kinetic system has multiple 
equilibria in certain ranges of the parameters that specify it. This is because of the following 
result concerning the "hydro dynamic" system (2.14). 

Proposition 2.3 There are choices of the parameters a, b, c and v for which there are 
nonempty intervals (Eq, E\) such that, if \E\ is outside of the closed interval [-Bo, E{\, then 
there is unique stable fixed point for the system (2.14), while, if \E\ e (Eq,Ei) then there 
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are three fixed points for the system (2.14). Moreover two of them are stable and one is 
unstable. 

Stability here is meant in the sense that the eigenvalues of the differential have a 
strictly negative real part. 

The proof of Proposition 2.3 is an explicit calculation which we omit (see [7] ) . However, 
to gain an intuitive understanding of why there should be multiple stable equilibria for 
certain parameter ranges, think of the e = limit of (2.1) as a constrained motion on the 
"manifold of Maxwellians" . Without this constraint, which is generated by the collision 
kernel, the evolution would be the one ruled by the electric field and the linear operator 
L. Clearly, this evolution has an unique attracting fixed point, which, for a > 0, is not 
Maxwellian. What happens is that there are one or more places on the non-linear constraint 
manifold that are locally closest to the attracting point of the unconstrained system. Each 
of these is a stable equilibrium for the constrained evolution. As the parameters are varied, 
the position of the unconstrained fixed point relative to the manifold of Maxwellians varies, 
and with this variation in geometry, the number of locally closest points varies. 

The main physical issues regarding this model are sttled at this point: We have proved 
the existence of the multiple stable equilibria at the kinetic level - e small, but positve 
- that had been found and investigated in [7] at e = 0. Moreover, we remind the reader 
that we do not know how to establish such a result using conventional expansion methods: 
the difficulty being that if E is not small, and hence possibly out of the range where 
multiple equilibria exist for (2.14), we only know how to prove (2.21) locally in time. This 
is insufficient to show that for e small enough, one never wanders far from any M( u * >e *) 
with (w*,e*) stable for (2.14). While the entropy methods we use do let us do this for 
arbitrary E, there are several finer question that one could ask, but for which we have only 
incomplete answers. 

First, one can ask whether or not there is an actual stationary solution inside the 
invariant neighborhoods that we have found of the M( u * jC »), and second, once one knows 
that stationary solutions exist, one can ask whether or not solutions actually tend to 
converge to one of these stationary solutions as t tends to infinity. 

The first question we can answer comletely only in the BGK case. The positive answer 
is given by: 

Theorem 2.4 

Let Q = Qbgk and a > 0. Then for each fixed point of (2.14) there is exactly one 
stationary solution of eq. (2.1). This solution lies in a suitably smalll neighborhood of 
M u *,e* , the Maxwellian corresponding to («*, e*), and it inherits the stability properties of 
the hydrodynamical fixed point. 

For the Boltzmann kernel we have only a partial result that is reported in Section 8. 
Also on the question of convergence we have only very partial results. These are reported 
in Section 9, where the difficulties are explained as well. But though it would be desirable 
to have a more complete resolution of these issues, they are not central to establishing 
that the kinetic systems does actually have the several stable regimes that one sees in the 
limiting "hydrodynamic" equations. 
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The proofs are organized as follows. In Section 3 we prove moment bounds. Section 
4 contains the proof of two "interpolation inequalities" . The first of these will be used to 
obtain a-priori smoothness bounds in Section 5. The second will be used to transform the 
smoothness bounds of Section 5 into a lower bound on the variance of our density. The 
smoothness requires Sobolev estimates for the collision kernel, which are straightforward 
for Qbgk, while for Qb they rely on some recent results [13], which we simply state here. 
Having assembled these moment, smoothness and interpolation bounds, we can use a key 
entropy production inequality for Qb proved in [11]. The analogous inequality for Qbgk 
is proven in a simpler way in Section 6. This is used to get quantitative bounds on the 
tendency of the collision operator to keep the density nearly Maxwellian. What we obtain 
directly is L 1 control on the difference between / and My, but the smoothness bounds 
together with the interpolation bounds allows us to obtain control in stronger norms. 
Section 7 contains the proofs of the Theorems 2.1 and 2.2, which, given the lemmas, are 
quite short. Section 8 is devoted to the proof of existence of stationary solutions, In section 
9 we discuss the tendancy toward these stationary solutions. 

3. Moment Bounds 

In this section we establish a-priori moment bounds for solutions of (2.1). In estimat- 
ing the evolution of the moments, we shall use one set of methods to treat the effects of 
collisions, and another set to treat everything else. Thus it is natural to rewrite (2.1) as 

-f( v ,t) = Cf(v,t) + -Q(F)(v,t) (3.1) 

where Cf = —E ■ V/ + Lf. 

We shall use the standard "bracket notation" notation for averages: ((f)) t denotes 
J" 3 <f)(v)f(v, t)d 3 v for any positive or integrable function </>. 

Throughout this paper, K will denote a computable constant that depends at most 
on the electric field E, the paramters a, 6, c and v specified in (2.4), and where indicated, 
also on the fourth moment of the initial distribution: (|t>| 4 )o. The constant will, however, 
change from line to line. 

In these terms, the main result of this section is: 

Theorem 3.1 Let f denote a solution to (2.1). Then for Q = Qb or Q = Qbgk, we have 

(\v\*) t <(\v\\ + K . 

Note that by Jensen's inequality, this imediately controls all lower order moments as 
well. The first step, however, is to directly conrol the second moments. 

Lemma 3.2 Let f denote a solution to (2.1). Then for both choices of the collision kernel 
Q, we have 

(\v\ 2 )t<(\v\ 2 ) + K (3.2) 
8 



Proof: Since (\v \ 2 )t is a collision invariant, and J v 2 L2fd 3 v = 0, 



> |2 >-L 



M 2 £,/(M)d 3 t> 

+ 2E ■ (v) t - 2 [ D(v)M(v)v • V (^4) d 3 v = 
Jr3 V. M(v) J 

2E ■ (v) t + 6(D) t + 2(v ■ VD)((v ■ VD)) t - 2(D\v\ 2 ) t 
Now observe that 

(v ■ V.D) < (3.3) 

and that 

c<D<(a + c) (3.4) 
for all v. Finally, by Jensen's inequality, 

l^> t |<«M 2 >t) 1/2 (3-5) 
These facts, combined with the previous calculation, yield the estimate 

^(M 2 >* < -2c(|t;| 2 ) t + 6(a + C )+2|i?|((|^| 2 ) t ) 1/2 
Straightforward estimation now leads to 

±(\v\ 2 ) t <-c(\v\ 2 ) t +K . 

Then (3.2) in turn follows from the fact that any solution of the differential inequality 

x(t) < -cx(t) + K satisfies x(t) < x(0) + K/c. □ 

We next parley these bounds into bounds on the fourth moments; i.e., (H 4 )t- Since 
(|w| 4 )t is not a collision invariant, these depend on the particular collision kernel under 
consideration. 

Lemma 3.3 Let Q = Qbgk ■ Then for any density f , 

[ \v\ 4 Q(f)d 3 v<(160/3)((\v\ 2 ) t ) 2 -(\v\% 

Proof: By an easy calculation, 

/ \v - (v) t \ 4 M f (v)d 3 v = 
9 



(5/3) (L \ v -^\ 2M f^ d3v ) = ( 5 /3)(k- (v)t\ 2 ) 2 t 

Next, note that \v\ < \v — {v) t \ + \{v) t \, and thus, \v | 4 < 16(\v — (v) t \ 4 + \(v)t\) 4 - Combining 
this and Jensen's inequality as in (3.5) with the above, we have 

/ \v\ 4 M f (v)d 3 v < (160/3) (<|v| 2 > £ ) 2 . 

The result now follows directly from the form of Q(f). □ 

Lemma 3.4 Let Q = Qb- Then there are positive constants c\ and C2 such that for any 
density f , 

I \v\ 4 Q(f)d 3 v < Cl ((\v\ 2 ) t ) 2 - c 2 (\v\% 
Jtr 3 



Proof: The result follows once we prove that 

J d 3 v\v\ 4 fof(v) < (£-c 2 )(\v\ 4 ) t + Cl ((\v\ 2 ) t ) 2 . 

To check this, let 7r w and tt^j denote the projection in the direction u> and the complemen- 
tary projection respectively. Then we can write v as 

and 

\v\ = \Tr u v*\ + |?r v„\ . 

Hence 

\v\ 4 = ItTuV*] 4 + |tt^X| 4 + 2|7r w ^| 2 |7r^-^| 2 . 

Averaging on S 2 we get 

/ duB(\v-v*\,u)[\n LJ v*\ 4 + |tt^<| 4 ] = |^| 4 tt / d6\ sin9\h{9) [cos 4 9 + sin 4 9] 

Js 2 JO 

= \v\ 4 (£-2tt J d9\ sin 9\h(9) cos 2 9 sin 2 9^j = \v\ 4 (£-c 2 ). 

□ 

Proof of Theorem 3.1: Calculating as before, we have: 
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d /I |4 



- f \v\ 4 Cf(v,t)d 3 v + - [ \v\ 4 Q(f)(v,t)d 3 v = 

4E-(\v\ 2 v} t -4 [ D(v)M(v)\v\ 2 v-v(^f]d 3 v + - [ \v\ 4 Q(f)d 3 v = 
in 3 V M ( v ) J e Jn 3 

4E ■ (\v\ 2 v) t + 12(\v\ 2 D) t + 4(\v\ 2 v ■ VD)((v ■ VD)) t - 2(D\v\ 4 ) t + - I \v\ 4 Q(f)d 3 v 

6 J-R 3 

We now estimate using (3.3) and (3.4) just as in the proof of Lemma 3.2, together 
with Lemmas 3.3 and 3.4, to control the collision term. (Note that once again L 2 makes 
no contribution.) The result is 



i\E\ ((M 4 ),) 3/4 + 12(o + c)(M 2 }i - 4 c (|t,| 4 >, + i(c, ((|f I 2 ),) 2 - c 2 (|f| 4 )i) • 
Next, by Lemma 3.1 together with Jensen's inequality, we know that 12(a + c)(\v \ 2 ) t is 

1 /2 

bounded above by a universal constant plus | 4 )o) • Thus, if we introduce the scaled 
time parameter 

r = (l/e)t , 

and define x(r) := (M 4 )t, we have that x(t) satisfies a differential inequality of the form 

x(t) < -x(r) +K 

from which the result follows. □ 
4. Interpolation inequalities 



The lemmas in this section are several interpolation inequalities related to the familiar 
Gagliardo-Nirenberg inequalites, but with some special features adapted to our applica- 
tions. The inequality of Lemma 4.2 is the most novel and interesting of these. 

Lemma 4.1 Let f e L^IR 3 ). Then there is a universal constant C such that if f has a 
square integrable distributional Laplacean, then f has a square integrable gradient, and 

||V/|| 2 <C||/||? /7 ||A/||5 /7 (4.1) 

Similarly, there is a universal constant C such that if A 2 f is square integrable, then Af 
is also square integrable, and 

IIA/lla^CH/llJ^IIAVll^ 11 (4-2) 
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Proof: Taking Fourier transforms, we have 



liv/ni 



\p\ 2 \f(p)\ 2 d 3 P 

|2|JV„m2j3„ 



P i<i? 

Ci2 5 ||/||? + iT 2 
Ci2 5 ||/||? + i2- 2 ||A/||l 



|pH/(p)| 2 d 3 p< 



P |>i? 
p| 2 |/(p)| 4 d 3 p < 



where the computable constant C changes from line to line. Optimizing over R now yields 
(4.1). The proof of (4.2) is done in the same way. □ 

The next inequality is similar in effect to an "uncertainty principle" . We shall use it 
to obtain uniform lower bounds on the variance of our density /. 

Lemma 4.2 Let f e L 1 (M 3 ), and suppose that f has a square integrable distributional 
gradient. Then there is a universal constant C such that 



R 3 



f(v)d*v < C||V/||* 



4/9 



R 3 



V — 



R 3 



vf(v )d s v 



f(v)d 3 v 



5/9 



(4.3) 



Proof: The right side of (4.3) is decreased when we replace / by its spherically symmetric 
decreasing rearrangement, while the left side is unchanged. We may therefore assume 
without loss of generality that / is spherically symmetric and radialy decreasing. Now fix 
R > 0, and define g(v) := f(v) for \v\ < R, and g(v) = otherwise. Define h by / = g + h. 
Now clearly, 

/ h(v)d 3 v < R- 2 [ \v\ 2 f(v)d 3 v . 

JtR 3 J\v\>R 

Also, 

/ g(v)d 3 v = 

JR 3 

(An/3)R 3 f(R) + (An/3)R 3 ((4n/3)- l R- 3 [ (f(v) 

\ J\v\<R 

(An/3)R 3 f(R) + (4n/3)R 3 ( (47r/3)" 1 J R- 3 f (f(v) 

\ J\v\<R 

Now, since / is monotone, 

/ |^| 2 /(^)d 3 ^ > f(R) [ \v\ 2 fd 3 v = f(R)(4n/b)R 5 . 

J\v\<R J\v\<R 
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Then, with A denoting the principle eigenvalue for the Dirichlet Laplacean in the unit ball, 
we have that 

/ {m-fiRrftfvKtfX- 1 [ \V(f(v) - f(R))\ 2 d 3 v < R'X-'WVfg . 

J\v\<R J\v\<R 

Combining the above, we have 

/ f(v)d 3 v<c(\\Vf\\ 2 R 5 / 2 + R- 2 [ \v\ 2 f(v)d 3 v) 
Jtr 3 \ Jn 3 J 

Optimizing over R yields the result. □ 

5. Smoothness bounds 

The purpose of this section is to establish a-priori smoothness bounds for solutions 
of (2.1). The main result is the the following: 

Theorem 5.1 Let f be a solution of (2.1) such that || V/(-, 0) || 2 is finite. Then, if Q = 
Qbgk, there is a constant K such that 

||V/(.,t)|| 2 <K(l+||V/(.,0)|| 2 ) (5.1) 

for all t > 0. Similarly, suppose that || A/(-, 0) 1 1 2 is finite. Then there is a constant K such 
that 

||A/(.,t)|| 2 <K(l + ||A/(.,0)|| 2 ) (5.2) 

for all t > 0. 

If Q = Qb the same results hold provided ||/(-, 0) — My(. j0 )||i and e are both sufficiently 
small. 

Remark: The smallness condition on ||/(-, 0) — Mj(. )||i poses no problem here since 
we are avoiding an initial layer by assuming that /(•, 0) = M^.^y However, it seems likely 
that it would be strightforward to include an inital layer analysis, and to dispense with this 
condition - even the present proof does not require ||/(-,0) — Mf(. j0 )||i to be particularly 
small. 

Proof: Once again, we write (2.1) in the form (3.1). Then differentiating, and integrating 
by parts, we have 

-2/ Af(v,t)Cf(v,t)d 3 v-2 [ Af(v,t)Q(f)(v,t)d 3 v . 

As with the moment bounds, we begin by estimating the individual contributions to 
the £ term. 
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Lemma 5.2 There is a universal constant K so that for all solutions f of (2.1), 

-2 [ Af{v,t)Cf{v,t)d 3 v < K 



for allt>0 

First, 2E ■ J M3 AfVfd 3 v = 0. 
Next, 

-2 [ Af(VD-Vf)d 3 v= [ |V/| 2 (A J D|d 3 t;<(a + C )||V/||| 

Now, repeatedly integrating by parts: 

-2 f Af(\7D-v)fd 3 v = 

2 / Vf(v(AD)f + 3VDf + (v ■ VD)Vf)d 3 v < 
Jr 3 

-3 [ f 2 (AD)d 3 v- [ f 2 (v-V(AD))d 3 v 
where we have once again used (3.3). This gives us a bound of the form 

-2 / Af(VD-v)fd 3 v<K\\f\\ 2 2 . 
We now use a standard interpolation inequality, the Nash inequality: 

<^liv/||^ /5 ||/||? /5 . 



This allows us to eliminate ||/||2 in favor of ||V/||2, the quantity of interest, and ||/||i, the 
conserved quantity. 

-2 / Af{VD.v)fd 3 v<K\\Vf\\T\\f\\t /5 ■ 
Jr 3 

-2 / Af(Dv ■ Vf)d 3 v = [ |V/| 2 V • (Dv)d 3 v < K\\Vf\\ 2 

JjR s JR 3 



Next: 



Apart from the collision term and the favorable dissipation term, the only other term to 
be bounded is —2 AfDf. However, integrating by parts, we obtain terms identical to 
terms that we have already bounded. 

The dissipation term is bounded using Lemma 4.1 as follows: 



-2 / D(Af) 2 d 3 v < -2c||A/|| 2 < 

^iiv/iij 4/5 ii/iir 4/5 = ^iiv/nj 4/5 
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since ||/||i = 1 for all times t. 

Thus, with x(t) := ||/(-, t) |||, we have established 

-2 / Af(v, t)Cf(v, t)d 3 v < K(-x(t) r/5 + K(x(t) 3/5 + x(tf /10 + x(t))) < K 

JTR 3 

since the largest power if x(t) has a negative coefficient. □ 

To conclude the proof of (5.1), we need smoothness bounds for the collision kernel; in 
particular, we need an estimate on the smoothness of the gain term in the collision kernel. 
For the BGK case, this is given by the following: 

Lemma 5.3: For any positive integer n there is a constant K depending only on the second 
moment of f , such that 

2 I (-Arf(v,t)Q(f)(v,t)d\<K-\\(-Ar/ 2 f\\t (5.3) 
Jr 3 



Proof: We only check the case n = 1; the rest are similar. We have 



(5.4) 



'R 3 

Now, since 



-2 / AfQ BGK (f)d*v = 
Jtr 3 

-2 / (AM / )/d 3 z;-2||V/||| . 



2 / (AMf)fd 3 v < CTjT 1 1| V/H2, 
Jn 3 



'R 3 

we can use Lemma 4.2 to get the bound 

TJ 1 < C\\Vf\\T (5.5) 

and hence 

-2 / (AM^/d^^C^HV/lb^CllV/II^^K+llV/lli. 

JR 3 

Combining this with (5.4) we get the Lemma 5.3 for n = 1 □ 

The analog of Lemma 5.3 for Qb is slightly more complicated. 

Lemma 5.4: For any positive integer n there are constants K and 5 depending only on 
the a-priori bound on the second moments of f , such that 

2 / (-Arf(v,t)Q B (f)(v,t)d s v < K - (l/2)\\(-Ay/ 2 f\\l (5.6) 

JlR 3 
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whenever 

||/(-,t)-M /( , t) (.)||i<5. 
where t is the constant in (2.13). 

Proof: It is sufficient to note that 

2 / (-ArfQ B (f)d 3 v = 
Jtr 3 

n [ f(-A) n f o fd 3 v - 2e\\(-Ar/ 2 ff 2 < 

Jtr 3 

2£(||(-A)"/ 2 /|| 2 ||(-A) n/2 / o f\\ 2 - ||(-A)"/ 2 /ll^ 

What we need now is control over ||(— A) n / 2 /o/|| 2 . This is provided by an inequality from 
[13], where it is shown that for any 7 > there are constant K and 6 > such that 

||(-A)-/ 2 / o /|| 2 < K + 7 ||(-A)-/ 2 /|| 2 . (5.7) 

whenver 

\\f-M f \\ 1 <5. 

This inequality is proved in [13] assuming that / has zero mean and unit variance, 
and under these conditions, the constants K and 5 are universal. Scaling the inequality, we 
have it holding with constants K and 6 depending only on the second moments of /. (The 
inequality is applied in [13] to get stong exponential convergence estimates for the spatially 
homogeneous Boltlzmann equation with physically realistic constants in the bounds.) 

With this inequality, we need only take 7 = 1/2. □ 

Proof of Theorem 5.1: We begin with the proof of (5.1). 

To put all of the lemmas together, let r := (l/e)t as before, and put x(r) := ||/(-, t) Hi- 
Combining Lemmas 5.2 and 5.3, we get in the BGK case that 

x(t) < eK + K — x(t) 

and any solution of this differential inequality satisfies x(t) < x(0) + (1 + e)K for all r. 
This establishes (5.1) in the BGK case. The proof of (5.2) in the BGK case is entirely 
analogous. 

To handle the Maxwellian collision kernel case, first define 

t = inf{t such that 7||/(-,t) -M /( . jt) || 2 > 5/2} (5.8) 

where 6 is the universal constant from Lemma 5.4. We take our initial condition small 
enough that i > 0. Then, for all t < i , we have the following differential ineqaulity by 
combining Lemmas 5.2 and 5.4: 

x(r) < eK + K — £x(r) 
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and any solution of this differential inequality satisfies x(r) < x(0) + (1 + e)K/£ for all 
r<(l/e)t. 

It only remains to show that actually t = +00. We shall do this in the next section 
using an entropy inequality. The entropy inequlity requires the a-priori smoothness bounds 
from Theorem 5.1, and shows that a long as they hold; i.e., as long as t < t, we have an 
upper bound of the form 

||/(-) t) — M/(. )t ) ||i < 

"entropy \t ) • 

All we have to do now is to take e so small that 

<5entropy(e) < 5/4 

and then it is clear from the definition of i that i = oo. Thus, borrowing the entropy 
bound from the next section, (5.1) is established for Q B for all e sufficiently small, and all 
initial data sufficently close to a Maxwellian. Again, (5.2) for Qb is handled in an entirely 
analogous way. □ 

6. Entropy bounds 

Let h(p\\p2) denote the relative entropy of two probability densities p\ and p2 on IR 3 : 

Here we are primarily interested in bounds on h(f\Mf). It will be convenient however, 
to first obtain bounds on h{f\M), and to then relate the two relative entropies. We do 
this in the next two lemmas. 

Lemma 6.1 Let f be any solution of (2.1) with 

(M 4 )o<C and ||V/|| 2 <C 
for all t > 0. Then there is a constant K depending only on C such that 

^-h(f(;t)\M)<K+- [ f(;t)Q(;t))d 3 V . (6.1) 



Proof: Differentiating, we have 

^h(f(;t)\M) 



[ %(lnf-lnM)d 3 v = 
Jtr 3 at 

[ (-E -V/ + Lf + e" 1 Q(/))(ln/ -lnM)d 3 v 
Jtr 3 



in 3 
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Now, integration by parts reveals that 

I Lf (In / -In M)d s v <0 
Jn 3 

and since InM is linear in \v\ 2 , 

- [ lnMQ(/) = 0. 

Finally, 

/ (E ■ V/) (in / — In M) d 3 v = 
Ju 3 

[ (E -Vf)lnfd 3 v - [ (E -Vf)lnMd 3 v = 
Ju 3 Jn 3 

E-(v) t 

Using Jensen's inequality to bound \(v)t\ in terms of the uniformly bounded (|f| 2 )t, we 
have the assreted result. □ 

Lemma 6.2 Let f be any solution of (2.1) with 

(M 4 )o<C and ||V/|| 2 <C 

for all t < T , some T > 0. Then there is a constant K depending only on C so that such 
that 

±h(f(;t)\M fi . tt) )<K+- f f(;t)Q(;t))d 3 V . (6.2) 

at e Jj^3 

for all t < T. 

Proof: By the definitions, we have 

h(f(.,t)\M f( ., t) )-h(f(.,t)\M) = 

(lnM) t - (lnM /( , t) )i = 

(3/2) ln(|^ - <^> t | 2 >, + (3/2) - (l/2)<|^| 2 > t 

Thus, 

±h(f(;t)\M f{ , t) ) = 

±h(f(;t)\M)+ 

(3/2)^\n(\v-(v) t \%-(l/2)^(\v\ 4 ) t 
Computing further with the logarithmic derivative term, 
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f f H\v - (v) t \ 2 ) t = 



To control this term, we need an upper bound on ((|t> — (v)t\ 2 )t) 1 But since (\v — 
(v)t\ 2 )t > C\\ V/H2 4 ^ 5 > K by Lemma 4.2, we have on application of Lemma 5.1 that 

^H\v-{v)^) t <K^({\v\ 2 ) t -{vf t ) 

The lemmas of Section 3 provide uniform bounds on this last term, and therefore provide 
uniform bounds on the derivative of h(f(-, t)\M^.^) — h(f(-, t)\M). The result now follows 
from the previous Lemma. □ 

Lemma 6.3 Let Q = Qbgk ■ For any density f with finite second moments, 

[ \nfQ(f)<-h(f\M f ) (6.3) 

Proof: Let S(f) = — J m3 flnfd 3 v denote the entropy of /. For any density / and for 
any r G [0, 1], put /(r) = (1 - r)f + rM f . We have 

±f(r) = M f -f = Q(f) 

and /(0) = /. Then 

/ lnfQ(f)d 3 v = -^-S(f(r)). 
Jtr 3 dr 

However, the entropy functional is concave, so that 

S(f(r))>(l-r)S(f)+rS(M f ) 

and hence 

s(/M)-s(/(Q)) s5(M/) _ g(/) 

Since Mf and / share the same hydrodynamic moments, S(Mf) — S(f) = h(f\Mf) □ 
Lemma 6.4 Let Q = Qbgk and f be any solution of (2.1) with 

(M 4 )o<C and ||V/|| 2 <C 
for all t Then there is a constant K depending only on C so that 

^h{f{^t)\M n . A )<K--h{f{^t)\M f{ .^ . (6.4) 
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Proof: This follows immediately upon combining the last three lemmas. □ 

In the case of the Boltzmann collision kernel, Lemma 6.3 is replaced by the following 
proposition proved in [11]: 

Proposition 6.5 For all C > 0, there is a positive function $c( r ) strictly increasing in 
r, such that for all densities f with 



f 



\v\ 4 f(v)d s v < C and ||V/|| 2 < C 



[ \nfQ B (f)<-$ c [h(f\M f )]. (6.5) 
Consequently, for Q — Qb, (6.4) is replaced by 

^h(f(-,t)\M f{ ., t) ) < K - h c [h(f(-,t)\M f( . jt) )} . (6.6) 

Now if /(-,0) = My(. j0 ), so that h(f(-, 0)|My(. jO ) = 0, then it is evident from this 
differential inequality that 

<S> c [h(f(;t)\M f( ., t) )}<eK 
for all t > 0. This together with Kullback's inequality 

\\f(;t)-M f{ . jt) f 1 <2h(f(;t)\M f{ . jt) ) (6.7) 

clearly implies that there is a function d entTOp y(e), decreasing to zero with e so that 

\\f(-,t) -M /( . )t) ||i < 5 entro py(e) 

for all t < i -the times for which we know that f(-,t) will satisfy the bounds in the 
hypothesis of Proposition 6.5. But as explained at the end of the proof of Theorem 5.1, 
this is enough to show that t = oo for all sufficiently small e. Thus, the entropy bound 
provides the information needed to complete the proof of Theorem 5.1 as claimed, and 
moreover, (5.8) hold globaly in time. 

7. Proof of Theorems 2.1 and 2.2 



Proof of Theorem 2.1 Let us consider first the case Q = Qbgk- Since Maxwellian 
initial data satisfies the hypotheses of Lemma 6.4, we have that (6.4) holds for the solution 
f(-,t) of (2.1) under consideration in Theorem 2.1. But then 

^(e^h(f(.,t)\M f( . !t) ))<Ke^ . 



20 



Since by hypothesis, h(f(-, 0)|Mf(. jO )) = 0, we have 



h(f(;t)\M fi ., t) )<eK 

for all t > 0. This together with Kullback's inequality (6.7) yields the first inequality 
asserted in Theorem 2.1, with Si(e) = e 1 / 2 . 

To obtain the second, note that since /(•, 0) is Maxwellian, there is a bound on || A/H2 
depending only on (\v\ )q. Inequality (5.2) of Lemma 5.1 now fives us a uniform bound 
on || A/(-, t) H2. Combining this with the interpolation inequality (4.1) finally yields the 
second inequality of Theroem 2.1 with 62(e) = e 1 / 7 for Q = Qbgk- 

The proof for Q = Qb is only slightly more involved, but in fact we have already given 
the proof of the first part of Theorem 2.1 in our "back and forth" proof of the smoothness 
bounds and entropy bounds for this case. As observed at the end of Section 6, 

\\f(-,t) -M /( . )t) ||i < 5 en tropy(e) 

for all t > 0, and 5 en tro P y(e) does decrease to with e as required. The second part follows 
in an entriely similar way. □ 

Proof of Theorem 2.2 This will follow from Theorem 2.1, and it is now no longer 
necessary to separate the cases Q = Qbgk and Q = Qb- 
As we have computed in the proof of Lemma 3.1, 



t — 

2\ 



- 2E ■ (v) t + 6(D) t + 2(v ■ VD)((v ■ VD)) t - 2(D\v\ 2 ) t 

If we replace the density /(•, t) everywhere on the right by Mj(. )t ), by definition we obtain 
the function G(u(t), e(t)) where u(t) and e(t) are the moments of /(•, t) figuring in Theorem 
2.2. The error we make has to be estimated term by term. The least trivial of these terms 
concerns the contribution from (D\v\ 2 ) t , and is estimated as follows: 

\(D\v\ 2 ) t - [ D(v)\v\ 2 M f( . :t) (v)d 3 v\ < 

[ D(v)\v\ 2 \f(v,t)-M fi ., t) (v)\d 3 v < 
/ D(v)\v\ 2 \f(v,t)-M f( ., t) (v)\d 3 v + 

J\v\<R 

[ D(v)\v\ 2 \f(v,t)-M f{ . jt) (v)\d 3 v < 

J\v\>R 

C^ 2 ||/(-,t)-M /( , t) || 1 + i?- 2 ^J^| 4 |/(z;,t)-M /( , t) (z;)|d 3 ^ < 
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C(i? 2 ||/(-,t)-M /( , t) || 1 + J R- 2 K((|z;| 2 ) t + l)) 
Optimizing over R now yields the result 

\(D\v\*) t - [ ^)|^| 2 M /( . jt) (z;)d 3 t;|<K||/(-,t)-M /( ., t) ||; /2 

Theorem 2.1 yields a bound of the size 51(e) 1 / 2 . The other error terms in the time deriva- 
tives of (H 2 )t and (v)t are bounded by a direct application of Theorem 2.1 (and hence 
yield errors of order Si(e) instead of ^i(e) 1 / 2 ). □ 

8. Stationary solutions 



In this section we discuss stationary solutions of (2.1), and prove Theorem 2.4. The 
stationary solutions of (2.1) are the positive normalized solutions of 

eCf + Q(f)=0. (8.1) 

The existence of such solutions relies on a simple fixed point argument. Recalling the 
expression of Q(f), we put 

Mf in the BGK case , 
/ o / in the Boltzmann case . 



J(f) = 

Equation (8.1) can be rewritten as 



/ = Hf) + jCf, (8.2) 



where £ = 1 in the BGK case. The explicit form of the linear operator C implies that, for 
e sufficiently small, the operator 

Q={l-\C)-\ (8.3) 

is positivity and normalization preserving on Li(IR 3 ). (That is, it is a Markovian operator). 
Moreover, it is clear from the smoothing properties of of Q that it is compact on Li(IR 3 ). 
Also, J is a positivity and normalization preserving map of Li(IR 3 ) into itself. Then since 
we can rewrite (8.2) as 

f = GJ{f). (8-4) 

we see that the solutions of (8.1) that we seek are the fixed points of the map / i— > QJ(f). 
The properties of this map listed above prove the existence of fixed points (see [14]). Let 
us denote by /* one of them. It is then easy to check that it solves (8.1) pointwise. 

Next we note that the same arguments used in previous sections imply that, if / 
solves (8.1) then 

||/* < 5 en tropy(e). 
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This concludes the first part of Theorem 2.4. 

Everything done so far would apply in the Boltzmann well as the BGK case. 

What we don't know at this point is: Are there fixed points in each of the stable neigh- 
borhoods, and if so is there exactly one in each stable neighborhood. 

These questions can be positively answered in the BGK case in a simplre way. It 
seems likely that one could also provide a positive answer for at least the first of them for 
the Boltzmann kernel, but we have not done more than skecth a lengthly argument, and 
so will confine ourselve to the BGK case. 

In fact, since Q = I + eCQ with / the identity map and J(f) — Mf, we can write (8.4) 

as 

f = M e + eCgM e . (8.5) 

Note that the right hand side of (8.5) depends only on u and e, the first and second 
moments of /. If we multiply (8.5) by v or by v 2 and integrate, we get 

F e (u,e) = 0, , . 

8.6 

G e (u,e) = 0, 

because / and Mf have the same first two moments. The functions F e and G e are quite 
complicated, but for e = they reduce to the functions F and G in the right hand side of 
(2.14). Then, by Proposition 2.3 we know there are solutions (it*,e*) to (8.6) for e = 0. 
Moreover the differential of the map (u, e) — > (F e (u, e), G(u, e)) has eigenvalues with non 
vanishing real part, in e = and (u,e) = (tt*,e*), when E is in the appropriate range. 
Therefore, by the implicit function theorem, for e sufficiently small, we have a unique 
solution (u e ,e e ) in a neighborough of (tt*,e*) to (8.6). Let M e be the Maxwellian with 
moments (u e , e e ). Then it is easy to check that 

/ = M e + CQM e 

is solution to (8.1). This concludes Theorem 2.4. 

9. Long time behavior. 



It is natural to ask whether the stationary solution are the asymptotic limits as t — > 
+oo of the evolution starting in appropriate neighborhoods of the fixed point Maxwellian. 
To this we have only a partial answer even in the BGK case: 

Proposition 9.1 Choose a stable fixed point (u*, e*) of (2.14) and let /* be a station- 
ary solution of (2.1) in the neighborhood of M^ u * e *y Assume that the solution ft of the 
time dependent problem starting near M( u * )6 *) has moments u(t) and e(t) converging to 
and e* respectively. Then 

lim ||/-M| 2 = 0. 

Unfortunately we do not have enough control on the time behavior of the solution 
to check the convergence of the moments. We expect however such convergence and this 
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can be proven for a modified model where we consider, instead of a diffusion coefficient 
D depending on the velocity, one depending only on the average e. A straightforward 
calculation then shows that one gets closed equations for the first two moments and the 
long time asymptotics is easily obtained. In this case all our results still apply and the 
conditions of Proposition 2.5 are fulfilled. 

First we prove Proposition 2.5. Let /* be a fixed point and / the distribution at time 
t. Calculating as in previous sections, we have: 



1 d 

2dt 



'R 3 

It is easy to chech that 



11/ - Mil = / d 3 v(f(v) - f*(v))[Cf(v) + -Q(f)(v)} 

JR 3 e 

d 3 v(f(v) - Mv))[C(f(v) - Mv)) + -(Q(f)(v) - Q(f*)(y))] 

R 3 e 

= / d\(f(v)-Mv))C(f(v)-Mv)) 
+ - [ dMf(v)-Mv))(M f (v)-M f Av))--\\f-M\l 



'R 3 

On the other hand 

d a v(f(v) ~ Mv))Lt(f(v) - f.(v)) 



d s v(f(v)-f*(v))L2(f(v) -f.(v)) <0, 

R 3 

/ d 3 v(f(v) - Mv))E ■ V v (f(v) - /*(«)) = 0. 
Jtr 3 



R 3 



d 3 W v (f(v) ~ U{v))D{v)M{v)vJ HV) 1 ! '° 

R 3 



M 



= - [ d 3 v\V v (f(v) ~ Mv))\ 2 D(v) - [ d 3 vv ■ V v (f(v) - Uv))D(v)(f(v) - /.(«)) 

iR 3 iR 3 

= - / d 3 v\V v (f(v) ~ f*(v))\ 2 D(v) - \ f d 3 vv ■ V v (f(v) - f*(v)) 2 D(v) 

JR 3 Z JR 3 

= - / d s v\V v (f(v) - f*(v))\ 2 D(v) + \ ! d 3 v(f(v) - Uv)) 2 V v ■ (vD(v)) 
Jr 3 1 Jr 3 

= - / d 3 v\V v (f(v) ~ h(v))\ 2 D(v) + 1 - f d 3 v(f(v) - U(v)fv ■ V v D(v)) 

JR 3 1 JR 3 

+ |/ d 3 v(f(v)-Mv)) 2 D(v) 



'R 3 

<l I d 3 v(f(v)-Mv)) 2 D(v) 

'R 3 



2 
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The last ineqality is consequence of the (3.3). Since 

f d 3 v(f(v) - f.(v))(M f (v) - M fm (v)) <\\\f-U\\l + \\\M f -M ft HI, 
we conclude that 

^11/ - fA\l < - e \\Mf - M u \\l ~ [\ ~ 3(a + c)]||/ - M\l 

With x(t) = ||/ — /* || I an d a = \\Mf — ||| we have, for e sufficiently small: 

x + cx < a 

and hence 



x(t) < xoe- 1 ^ + [ dse- {t - s ^ ea{s \ 
Jo 



This implies x(t) — > as t — > oo, provided that a(t) — > as t — > oo. This concludes 
Proposition 2.5. 

The convergence of My to is not easy to get for (2.1). A simple answer is obtained 
if one replaces the operator L\ given by (2.3) with 

LJiv) = D(e)V ■ (^)V/(-^)), (9-1) 

This model is much simpler than the one already considered, but still has a non trivial 
behavior on the hydro dynamical sale. In particular, it is easy to see that, if we write the 
equations for u and e, we get a closed system in those variables. Namely, (2. 14), (2. 15) and 
(2.16) are replaced by 

d( u (t)\_(F(u(t),e(t))\ 

the functions F and G being given explicitly by 

F(u,e) = E-u[v + D(e)\, (9.3) 

G(u,e) = Eu + 3D(e) -2eD(e). (9.4) 

Equation (9.2) is exact for this model independently of e and for suitable choices of the 
functions D(e) has several critical points for E in an appropriate range. Moreover the 
asymptotic behavior for large times is easy to establish. All our results apply to this 
model without substantial changes. In particular, in this case we can use Proposition 2.5 
to obtain the convergence to stationary solutions for large times. 
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